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Abstract 

We derive an effective ID theory from the Hamiltonian of the 3D system which consists 
of a mesoscopic conductor and reservoirs. We assume that the many-body interaction have 
the same magnitude in the conductor as that in the reservoirs, in contrast to the previous 
theories which made the ad hoc assumption that the many-body interaction were absent in 
the reservoirs. We show the following: (i) The effective potentials of impurities and two- 
body interaction for the ID modes become weaker as x goes away from the conductor, (ii) 
On the other hand, the interaction between the ID and the reservoir modes is important in 
the reservoir regions, where the reservoir modes excite and attenuate the ID modes through 
the interaction, (hi) As a result, the current I\ of the ID modes is not conserved, whereas 
the total current I is of course conserved, (iv) For any steady state the total current 
(/), its equilibrium fluctuation (5I 2 ) eq at low frequency, and non-equilibrium fluctuation 
{5I 2 ) noneq at low frequency, of the original system are independent of x, whereas (5I 2 ) eq 
and (SI 2 ) noneq at higher frequencies may depend on x. (v) Utilizing this property, we can 
evaluate (/), (5I 2 ) eq , and ^I 2 ) noneq at low frequency from those of the ID current I\. (vi) 
In general, the transmittance T in the Landauer formula should be evaluated from a single- 
body Hamiltonian which includes a Hartree potential created by the density deformation 
which is caused by the external bias. 
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1. Identification of the single-body part of the 3D Hamiltonian 

We consider a system which consists of a mesoscopic conductor and reservoirs (Fig. 
1), in which electrons are interacting each other according to the three-dimensional (3D) 
Hamiltonian (in the Heisenberg picture), 



H = m + v Q 



d 3 r H° (r,t) + ^ 



1 

2 

h 2 



2m 



d 3 r J d 3 r' V (r,r',t), 
V 2 + u c (r) +u i (r) +e(j) ext (r) 



V(r,*), 



Vo(r,r',t) = p(r,t)v(r-r')p(r',t), 
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where ip(r, t) is the 3D electron field, u c the confining potential, u l the potential of impurities 
and/or defects, (f> ext the external electrostatic potential, and v(r — r') = e 2 /|r — r'| denotes the 
Coulomb potential between the local charge density p(r,t) = i/j<(r,t)i/)(r,t) — pbg(j)- Here, 
e0t(r, t)V'(r, t) is the charge density of the electrons and —epsc is that of the background 
charges which are assumed to be fixed. Note that we cannot simply add to H the position- 
dependent "chemical potential term" like — / d 3 r p(r)i/j^(r,t)i/j(r,t), because it does not 
commute with H. 

We are interested in a steady state, which may be either equilibrium or non-equilibrium. 
Its density operator is denoted by (, which is to be determined self-consistently. The average 
charge density of this state is Tr (ep(r, t) = ep av (r), which is a function of the bias voltage. 
Although p av (r) vanishes at r — > oo, it is finite and r dependent in finite regions of space, 
because mesoscopic systems (conductor plus reservoirs) are spatially inhomogeneous. The 
finite p av induces a long-range effect, which does not vanish by the screening, through the 
Coulomb interaction. To construct a meaningful ID theory, in which the interaction of a 
"ID field" (Eq. (fPTD) with a "reservoir field" (Eq. ([16])) is not strong, we must remove the 
long-range effect from the interaction Hamiltonian. This can be accomplished by rewriting 
Eqs. (D-@) as 

H° + V + V av 
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d 3 r H°(r,t) + 



d 6 r I d 3 r' 



V(r,r',t) + V a 
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-— + u c (r) + u\r) + &, 
2m 



V(r,r',t) = 6p(r,t)v(r-r')6p(r',t), 
V av (r, r') = -p av (r)v(r - r')p av (r'). 



Here, 5p(r, t) = p(r, t) — p a v(r) denotes the charge fluctuation around p av , and 
average electrostatic potential; 



J d 3 r'v(r - r')p av (r') + ec. 



ext 



r . 



(6) 

(7) 

(8) 
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>) av {r) is the 
(10) 



The advantage of rewriting H in the form of Eqs. (|5|)-(|rO"D is that effects of V is much 
weaker than those of the original interaction Vq. In particular, V does not cause the long- 
range effect because it is already included in H°. It is therefore much better to start with 
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these equations, taking H° (rather than Hq) as the single-body part. For example, we can 
derive the Landauer formula by neglecting V in Eq. (|5]); 

G=(2e 2 /h)T[(f> av }. (11) 

Here, the transmittance T[4> av ] is calculated from 7i°, Eq. ([/]), which includes the self- 
consistent potential 4> av . On the other hand, if we neglected Vq in Eq. ([]]), we would obtain 
G = (2e 2 /h)T[(p av = 0]. This is different from Eq. (|ITD when the bias voltage is finite. Since 
actual measurements on mesoscopic conductors and on the quantum Hall effect are often 
performed under a non-negligible bias, one should use Eq. (|ll|) for these experiments. 

2. Decomposition of the electron field 

Using the above form of the Hamiltonian, we now decompose the electron field into two 
parts. We are interested in the case where the variation of W(x) of Fig.l is slow, i.e., 
W(x) ~ W(x + Xp), where is a length of the order of the Fermi wavelength. Otherwise, 
undesirable reflections would occur at the boundary regions between the conductor and 
reservoirs, and experimental results for such samples would not be very interesting. We 
do not want reflections either in the conductor except for reflections by impurities and/or 
defects. If the bias voltage is not too high, these conditions can be stated as the assumption 
that the x dependence of w c (r) + e<j) av (r) = u(r) is weak, i.e., u(x,y,z) ~ u(x + X F ,y,z). 
We also assume that only the lowest subband (see Eq. ( |i~4"D below) is occupied by electrons. 

Under these assumptions we can find approximate solutions of the single-body 
Schrodinger equation 

<p(r) = e<p(r), (12) 
for e not far from the Fermi energy, in the following form: 

ip^{y,z;x), (13) 

where C is the normalization length in the x direction, and ^(y, z; x) is the normalized 
eigenfunction belonging to the lowest eigenvalue u _L (x) of the following eigenvalue equation 
(in which x is regarded as a parameter): 

ip^iy, z; x) = u ± (x)(p ± (y, z; x). (14) 

The x-dependent wavenumber Kj.(x) is defined by Sk = h 2 Kk(x) 2 /2m + m ± (x), where 
denotes the eigenenergy of the state <£>fc(r). The Sk can be expressed as Sk = h 2 k 2 /2m + e ± , 
where k = Kk(0) and e 1 - = m _l (0) are the wavenumber and the lowest subband energy, 
respectively, at the center of the conductor. We call y^'s "conductor modes". There are 
other solutions to Eq. (0), which are denoted by (p v . They include states whose e is far from 
the Fermi energy, states belonging to higher subbands, and states which decay exponentially 
in the conductor. We call these solutions "reservoir modes." The r dependence of the electron 
field ip(r, t) can be expanded in terms of ^'s and cp^s. We can thus decompose the electron 
field into two parts as ^(r, t) = i>c( r , t) + ipn(r, t), where 
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Mv, t) = £ Mr) I dh'Mr'W, *) = £ ^(r)c fc (t), (15) 
^(r, t) = £ y^(r) | d 3 rV:(r>(r', t) = £ ^(r)d„(t). (16) 



We call %pn the "reservoir field." 

We now construct a ID field ipi(x,t) from the 3D field ipc{ r ,t) as 

^i(x,t) = J J dydz ip ± *(y,z;x) ^c(r,f) ~£-^=exp K k (x)d, 



c k {t). (17) 



This ID field has the same contents as ^7, because the inverse transformation can be 
accomplished as ipciy.t) = <p (y, z; x)ipi(x, t). In terms of ipi and ipR, the Hamiltonian can 
be recast as 

H = H° c + V c + V C r + K + % + V av (18) 

Here, Hq {H r ) denotes the single-body part of the ID (reservoir) field, Vc (Vr) is its 
mutual interaction, and Vqr denotes the interaction between the ID and reservoir fields. In 
particular, H c = J dx 7i c (x,t) and Vc = \ J dx J dx'Vc(x,x',t), where, 



H° c (x,t) ^ViM) 



h 2 d 2 



+ u (x) + u\(x) 



Mx,t), (19) 



2m dx 2 

Vc(x,x',t) = 5pi(x,t)v 1 (x,x')5pi(x',t). (20) 

Here, Spi(x, t) = tpl(x, t)ij)i(x, t) — (tpl(x, t)ipi(x, t)) is the density fluctuation of the ID field, 
and 



u\(x) = J J dydz \^{y, Z] x)\ 2 u\v), (21) 
Vl (x,x') = J J dydz J J dy'dz' \<p L (y, z; x)\ 2 v(r - r') \<p L (y', z';x')\ 2 . (22) 

The 3D impurity potential w*(r) is spatially averaged, with the weight ly? -1 "! 2 , to give the 
effective ID impurity potential u\(x). In the conductor region, u\(x) has the same order 
of magnitude as w*(r). In the reservoir regions, on the other hand, u\(x) becomes smaller 
because ip 1 - extends over the width W{x) in the y direction, and u\(x) is the averaged value 
of the random variable w*(r) over this wide area. As x — > ±oo, in particular, u\(x) — > if 
W(x) — > oo. 

We can say the same for the two-body interaction potential vi(x,x'). In the conductor 
region, the major effect of the integrals over y, z, y', z' of Eq. (|22"D is simply to smear out the 
singularity of v(r — r') at r = r' by averaging over the conductor width. In the reservoir 
regions, on the other hand, the interaction becomes extremely weak as we show now. When 
x or x' is in a reservoir region(s), Vx(x, x') is screened by the reservoir modes through Vcr- 
It is therefore better to rewrite Vc + Vcr in Eq. (|L8"D as Vq c + V cr , where Vq c is the screened, 
short-range interaction and V CR = Vcr + Vc — V c c - For large \x — x'\ it is clear that V c c 
becomes vanishingly weak as \x — x'\ is increased. For x ~ x', on the other hand, the 
weakness of V c c is never trivial because the screening is ineffective at a short distance. To 
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investigate this case, we consider the behavior of the original, long-range potential Vi(x,x'). 
It is easy to show that Vi(x, x') oc l/W(x) as x ~ x' — > ±00. Hence V\ — > if W(x) — > 00. 
Since the interaction potential of Vq c decays more quickly than v\, it also approaches zero, 
more quickly, as x ~ x' — > ±00. 

We have thus shown that the potentials of both the impurities and two-body interaction 
for the ID field vanish as x and/or x' go away from the conductor. (It should be stressed that 
W(x) — > 00 is necessary for this property; otherwise, the impurity potential and the short 
range part of the interaction would remain finite in all regions.) This may partly justify the 
ID models [2] in which a ID field is assumed to become free as x, x' — > ±00. However, note 
that V' CR will cause additional effects on the ID field: the reserver modes excite, attenuate, 
and renormalize the ID field [1]. 

3. Calculation of the total current from the ID current 

We now turn to observables. What is measured in most experiments is the current / 
which is given by I(x, t) = J J dydzJ x (r, t), where J x denotes the x component of the current 
density. From the continuity equation for J and p, we have 

_/ (X)t) + _A (X)t)=0 , (23) 



where A is the ID density of electrical charge; A(x, t) = J J dydz ep(r, t). For a steady state, 

£<iOM)> = -g<' 



either equilibrium or non-equilibrium, Eq. (E3) yields -4-{I(x,t)) = — 4:{A(x,i)) = 0. That 



is, the average current {I(x,t)) is independent of x. 

Note, however, that this is not the case for the current fluctuations. For example, consider 
the symmetrized correlation function of 5I(x,t) = I(x,t) — (I(x,t)), 

C H (x, t-t') = (5I(x, t)5I(x, t') + 5I(x, t')5I(x, t))/2, (24) 

which may be taken as a definition of the current fluctuation. Equation fl2"3] ) and the time- 
translational invariance yield 

— C n (x, t) = — [C IA (x, t) + Cai(x, -t) - C IA (x, -t) - C AI (x, t)} , (25) 

where we have introduced the (non-symmetrized) density-current correlation, Ci\(x, t—t') = 
(SI(x, t)A(x, t')) and C A i{x,t — f) = (A(x,t)SI(x,t')) . It is seen that the current fluctuation 
generally depends on the position x. 

However, we can show that the low-frequency component of its spectral intensity, 
(SI 2 )^ = f^° 00 e luJt Cu(x,t)dt, is independent of x. To show this, we note that any real sys- 
tems which possess the thermodynamic stability must have the "mixing property" ; namely, 
any correlation functions decay as \t — 1'\ — > 00. This property might be lost if one took an 
oversimplified model, such as integrable models, for H. However, our H contains both the 
conductor and reservoir modes as well as their interactions. It is quite reasonable to assume 
that such a complicated 3D model has the mixing property. Therefore, Cn(t) vanishes as t 
exceeds a "correlation time" r c . Then, for w < l/rc we find from Eq. (|25|) 

— (5I 2 )^ ~ ^(^ 2 )-=o = 2 [C IA (x, 00) + C AI (x, -00) - C IA (x, -00) - C AI (x, 00)] . (26) 
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The right-hand side vanishes because of the mixing property. We have thus obtained the 
theorem: For a steady state, the average current is independent of x, and the low-frequency 
(uj <C 1/t c ) component of the spectral intensity of current fluctuation is also independent of 
x, whereas higher-frequency components may depend on x. This theorem is quite general: 
It applies to any steady states, either equilibrium or nonequilibrium, of real systems and of 
any theoretical models which are realistic enough so that they have the mixing property. 

The theorem is very useful for our effective ID theory. Let us define the current of the 
ID field by 



$0m) j^J^i(z,m + h.c. 



(27) 



Note that this ID current is not conserved because of the interaction Vcr, which transforms 
the ID current into the reservoir current and vice versa. At first sight, this might seem to 
cause difficulties. Fortunately, however, the above theorem guarantees that we can evaluate 
(I) and (5/ 2 ) w ~o & t arbitrary x. If we take x = (the center of the conductor), I(x,t) ~ 
Ii(x,t) because the reservoir field does not contribute there. Therefore, we can evaluate 
(I) and (SI 2 )^^ by evaluating (Ii) and ) w ~o &t x — 0. This may partly justify the 
Landauer-like approach [3], in which one evaluates the DC conductance and nonequilibrium 
current fluctuation by evaluating (7i) and (5If ) w ~o using an effective ID model. 

We have thus shown that it is sufficient to evaluate the ID current I\. Therefore, we 
can project out the reservoir degrees of freedom, which should behave as the usual two- or 
three-dimensional electrons because the reservoirs are large, in the standard manner [4]. By 
doing so, we can derive either the generalized Langevin equation [4] for I x , or the Liouville 
equation for the reduced density operator [4], ( c = Tr R [£], where the trace is taken over the 
Hilbert space of the reservor field. The original 3D model of the conductor plus reservoirs 
has thus been reduced to the ID model. Detailed discussions and applications to various 
problems will be described elsewhere [1]. 
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Figure caption 



Fig.l The mesoscopic conductor and reservoirs. 
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Fig.1 A. Shimizu and T. Miyadera 



